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ABSTRACT 

The presence of dust in quasar absorbers, such as damped Lyman alpha 
(DLA) systems, may cause the background QSO to appear reddened. We inves- 
tigate the extent of this potential reddening by comparing the optical-to-infrared 
(IR) colors of QSOs with and without intervening absorbers. Our QSO sample is 
based on the Complete Optical and Radio Absorption Line System (CORALS) 
survey of Ellison et al (2001). The CORALS dataset consists of 66 radio-selected 
QSOs at z cm > 2.2 with complete optical identifications. We have obtained 
near-simultaneous B and K band magnitudes for subset of the CORALS sample 
and supplemented our observations with further measurements published in the 
literature. In total, we have B— K colors for 42/66 QSOs of which 14 have in- 
tervening DLAs. To account for redshift-related color changes, the B— K colors 
are normalized using the Sloan Digital Sky Survey (SDSS) QSO composite. The 
mean normalized B— K color of the DLA sub-sample is +0.12, whereas the mean 
for the no-DLA sample is —0.10; both distributions have RMS scatters ~ 0.5. 
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Neither a student's T-test nor a KS test indicate that there is any significant 
difference between the two color distributions. Based on simulations which red- 
den the colors of QSOs with intervening DLAs, we determine a reddening limit 
which corresponds to E(B — V) < 0.04 (SMC-like extinction) at 99% confidence 
(3<r), assuming that E(B — V) is the same for all DLAs. Finally, we do not find 
any general correlation between absorber properties (such as [Fe/Zn] or neutral 
hydrogen column density) and B— K color. The two reddest QSOs with DLAs 
in our sample have HI column densities which differ from each other by an order 
of magnitude and moderate gas-to-dust ratios as inferred from chemical abun- 
dances. One of these two QSOs shows evidence for strong associated absorption 
from X-ray observations, an alternative explanation for its very red color. We 
conclude that the presence of intervening galaxies causes a minimal reddening of 
the background QSO. 

Subject headings: ISM:general, galaxies:high-redshift, quasars: absorption lines, 
dust, extinction 



1. Introduction 
1.1. Red Quasars 

Although quasars are usually regarded as objects which are characterised by an ultra- 
violet (UV) excess, the existence of red QSOs has been known for more than two decades 
(e.g. Smith & Spinrad 1980; Bregman et al 1981; Stein & Sitko 1984). Despite their long 
history, red QSOs have continued to attract much interest in the literature, most notably the 
on-going debate to identify the primary cause of their extreme optical-infrared (IR) colors. 
Probably the most popular explanation for red QSOs is extinction due to dust at the systemic 
redshift, a conclusion that is gaining substantial support from the current generation of large 
QSO surveys (e.g. Richards et al. 2003; Hopkins et al. 2004; Glikman et al 2004; White 
et al. 2003). There are also convincing examples of individual QSOs which appear to have 
been highly reddened due to local (to the QSO) dust (e.g. Rawlings et al. 1995; Gregg et al. 
2000; Stern et al. 2003). If dust is indeed a common ingredient in QSOs, the corollary is that 
samples based on optically selected sources could potentially miss the reddened population, 
whereas radio-selected quasars would represent a more complete sample (e.g. Webster et al 



1 These data were obtained from the NTT on La Silla (ESO program 072.A-0014(A, B, C, D)) 



- 3- 



1995; Barkhouse & Hall 2001). Indeed, Webster et al. (1995) have claimed that up to 80% 
of QSOs may be hidden from the eyes of optical surveys due to dust extinction. 

However, some observations do not support the hypothesis of dust local to the QSO 
as the cause of reddening. For example, millimetre observations (Drinkwater, Wiklind & 
Combes 1996) of 11 of the reddest Webster et al. (1995) sources failed to detect CO which 
argues against the presence of host galaxy dust in these cases. Other possible sources of QSO 
'reddening' include the effect of host galaxy starlight and synchrotron radiation (Francis, 
Whiting & Webster 2000; Whiting, Webster & Francis, 2001; Benn et al. 1998). 

1.2. Dust Reddening Associated with Intervening Galaxies 

One factor that could have a significant impact on reddening QSO light that remains 
relatively unexplored is the presence of intervening dust located in foreground galaxies. The 
identification of a handful of red lensed QSOs lends credibility to this scenario, in which 
dust in the lensing galaxy reddens the background source (e.g. Gregg et al. 2000; Falco et 
al. 1999, Wucknitz et al 2003). The issue of dusty intervening galaxies is of interest not 
only in order to assess its contribution to the red QSO population, but is also of crucial 
importance for the study of quasar absorption lines. Quasar absorbers, such as the damped 
Lyman alpha (DLA) systems which correspond to the signatures of high redshift galaxies 
imprinted on the continuum of a background QSO, represent one of our best tools for studying 
the early galaxy population. Previous surveys for DLAs have almost universally relied on 
large optical databases of QSOs to find intervening absorbers. Since only relatively bright, 
optically selected QSOs have been used to identify DLAs, there has been widespread concern 
that dusty intervening galaxies may obscure the background QSO so that it drops below the 
survey magnitude limit. Such a selection effect would result in a biased view of high redshift 
galaxies. This general concern has been quantified in a number of theoretical works which 
have attempted to estimate the amount of reddening that a QSO may suffer and consequently 
what fraction of the DLA population, and their contribution towards the cosmic gas and 
metals budget, may have been overlooked (e.g. Fall & Pei 1993; Pei & Fall 1995; Pei, Fall 
& Hauser 1998; Churches et al. 2004; Vladilo & Peroux, 2005). However, until recently 
there had been little observational progress in quantifying the number of missed DLAs due 
to extinction effects. 

The first systematic program to quantify dust bias was undertaken by Ellison et al. 
(2001) by searching for DLAs in an optically complete radio-selected sample of QSOs. The 
Complete Optical and Radio Absorption Line System (CORALS) survey was based on a ho- 
mogeneous sample of 66 radio-selected, flat spectrum QSOs taken from the Parkes catalogue, 
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complete down to 0.25 Jy. Extensive optical follow-up observations resulted in identifica- 
tions for all sources. Therefore, subsequent spectroscopic searches for DLAs were based on 
a complete sample, ideal for determining the extent to which intervening galaxies obscure 
background QSOs. Ellison et al. (2001) determined that the neutral gas mass density of 
DLAs (I^dla) may have been previously under-estimated by at most a factor of ~ 2. How- 
ever, this work also found tentative evidence that Qbla and n(z) (the DLA number density) 
may be under-estimated for bright QSO samples with magnitude limits B < 19.5. Ellison 
et al. (2004) re-visited this issue and suggested that this could be understood in terms of 
an Eddington bias if the optical luminosity function deviates from a simple power law. To 
extend the original CORALS sample to lower redshifts, Ellison et al. (2004) studied Mg II 
absorption systems in an analogous radio-selected QSO sample. The main conclusion of 
that work was that the number density of intervening systems selected via Mg II absorption 
in the range 0.6 < z < 1.7 is in excellent agreement with magnitude limited samples (e.g. 
Nestor, Turnshek & Rao 2005). Most recently, Akerman et al. (2005) have measured the 
metal abundances of z > 1.8 CORALS DLAs and found that the weighted mean metallicity 
based on zinc (an undepleted proxy for Fe) is only marginally higher (but consistent within 
the error bars) than previous samples. 

Although the results from CORALS indicate that magnitude limited surveys (at least 
those complete to B ~ 19) have not missed a significant fraction of the gas or metals budget, 
this is not to say that DLAs are dust-free. In fact, it has been known for over a decade, 
based on the ratios of refractory and non-refractory elements, that dust is present in DLAs 
(Meyer & Roth 1990; Meyer & York 1992; Pettini et al. 1994, 1997). Understanding the 
level of extinction caused by this dust has so far yielded conflicting results. For example, 
Fall, Pei & McMahon (1989) and Pei, Fall & Bechtold (1991) found a significant difference 
between the spectral indices of QSOs with intervening DLAs and a control sample. However, 
this result has not been reproduced with larger, homogeneous QSO samples. Currently, the 
most stringent limit for reddening due to intervening DLAs has been obtained from fitting 
spectral indices to ~ 1450 QSOs from the Sloan Digital Sky Survey (SDSS): E(B - V) < 
0.02 for SMC-like extinction (Murphy & Liske 2004). Despite this sensitive limit, studies of 
spectral indices suffer from two drawbacks: they utilize limited wavelength coverage in the 
optical regime and have not yet been conducted on unbiased samples (i.e. they were taken 
from magnitude limited surveys). 
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1.3. Optical-Infrared Colors and Extinction Curves 

In this work, we take a different approach to studying the effect of dust extinction from 
intervening DLAs by using optical-IR QSO colors to search for a reddening signature. This 
technique has the advantage of using a wide wavelength baseline, from the observed frame 
B to K bands, over which to detect the effect of extinction. Moreover, by studying QSOs 
from the CORALS survey, we have a sample which is optically complete and in which all of 
the intervening high redshift DLAs have already been identified and metallicities measured. 

The observed reddening obviously depends largely on the assumed extinction law, al- 
though the QSO and galaxy (i.e. dust) redshift also play a part (as we demonstrate below). 
Despite its ubiquity in astrophysical environments, the properties of dust remain relatively 
poorly understood (see reviews by Mathis 1990 and Draine 2003) and reliable empirical 
extinction data are available for only the Milky Way (MW) and Magellanic Clouds. In par- 
ticular, data are sparse in the far-ultraviolet (FUV, Hutchings & Giasson 2001) and below 
the Lyman limit there is simply no flux with which to work. FUV wavelengths are particu- 
larly important for the present work, which considers the optical-IR colors of high redshift 
QSOs. The most often used parameterization of Galactic extinction is that of Cardelli, 
Clayton & Mathis (1989; hereafter CCM89). The advantage of the CCM89 formulae is that 
£(A) = A x / A v can be calculated for any value of Ry 2 . Another popular parameterization of 
Galactic extinction is that of Pei (1992), although this is for a fixed value of R v = 3.08. A 
comparison of these two extinction law fits is shown in Figure 1, with an Ry = 3.08 for the 
CCM89 law in order to match the Pei (1992) curve. It can be seen that at IR and optical 
wavelengths, the two curves are in close agreement. However, in the FUV (A < 1200 A) 
the curves diverge; this is purely due to extrapolation of the different fitting functions into 
a regime where no data were available. Whereas an extrapolation of the Pei (1992) parame- 
terization leads to a flattening of £(A), the CCM89 curve rises quickly at FUV wavelengths. 
In Figure 1 we also show the extinction curves of the SMC and LMC as fitted by Pei (1992), 
again extrapolated to the FUV using the published fitting functions. 

We illustrate the expected reddening due to intervening dust in Figures 2, 3 and 4. In 
each figure, we plot the theoretical reddening (as measured by a change in the B— K color) of 
a QSO at function of E(B — V). Since the color change depends on the redshift 

of the intervening dust, four different lines are plotted for dust at z& nst = 

2.0,2.3,2.6,2.9. The 

reddening induced by the two Galactic extinction curves is quite similar, although a CCM89 
law produces slightly redder colors due to the rapid rise of {(A) in the FUV. Reddening is 
most severe for SMC type dust, and could be even more extreme if the true extinction law 



2 In the usual notation, Ry = Ay/E(B — V), i.e. the total-to-selective extinction. 
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continues to rise at FUV wavelengths. That is, the parameterization of SMC dust by Pei 
(1992) may be considered conservative. 

2. Observations and Data Reduction 

Since flat spectrum radio-loud QSOs (RLQs) are highly variable, it is important to 
obtain optical and IR data from near-simultaneous epochs (e.g. Francis, Whiting & Web- 
ster 2000, hereafter FWW00). We achieved this by using the optical (SuSI2) and IR (SofI) 
imagers on the New Technology Telescope (NTT) at La Silla observatory. Since both in- 
struments are mounted on the same Nasmyth platform, the overhead incurred by switching 
between them is typically only about 2 minutes, depending on the source position relative 
to the zenith. Our strategy was therefore to observe a sequence of QSOs in the optical and 
then change to the IR. Optical and IR photometric standard stars were observed throughout 
the night at airmasses up to 1.8. 

Our data were obtained during two runs in October 2003 and March 2004. The March 
2004 run experienced excellent conditions: photometric transparency and typical seeing ~ 
0.6 arcsec. The earlier October 2003 run, suffered from somewhat poorer seeing (~ 1 arcsec) 
and the end of the last night (during which only SofI was used) was not photometric. 

2.1. Optical Imaging 

In addition to B-band images, we obtained supplemental V-band images for most QSOs 
with redshifts z cm > 3.0. Typical exposure times in the B and V bands ranged from 10 to 
200 seconds (see Table 1), the total integration was divided between two exposures with a 
spatial offset of a few hundred pixels between exposures. We used SuSI2 in 2x2 binning 
mode which yielded a pixel scale of 0.16 arcsec/pixel. Standard IRAF procedures were used 
for the reduction of the optical data: a bias was constructed from a median combination of 
five frames and the overscan was fitted and subtracted interactively. The initial flat-fielding 
was executed using dome-flats to correct pixel-to-pixel variations, followed by an illumi- 
nation correction of large-scale variations constructed from a super-sky flat. Instrumental 
magnitudes were determined for each night using IRAF's qphot package with between 20 
and 40 standard stars taken from Landolt (1992). The photometric solution (zero points 
and extinction coefficients) for each night in the B and V bands were determined by a least 
squares fit to the data. The zero points determined in this way indicate that conditions 
were photometric throughout the optical observations. The RMS error in the photometric 



-7- 



solution was 0.04 magnitudes in the V and B bands for the October 2003 run and 0.04 and 
0.07 for V and B bands respectively in March 2004. The final quoted B and V-band errors 
in Table 2 include the uncertainty in the fit of the photometric solution and photon statistics 
added in quadrature. 

2.2. Infrared Imaging 

Each SofI K s band integration consisted of a series of co-added 10 second exposures 
with a random offset within a 40 arcsec jitter box between each exposure. Total on-source 
integration times ranged between 100 and 5400 seconds (see Table 1). IR data reduction was 
carried out using Pat Hall's Infrared Imaging Reduction Software (PHIIRS) which incorpo- 
rates many standard IRAF routines. We did not follow the usual procedure of obtaining 
dark frames due to the non-linearity of dark current with incident flux. Instead, we used 
the PHIIRS task irsky to combine individual dithered exposures to simultaneously remove 
dark current and sky. We used 4-6 frames to construct the sky image for each science 
exposure; fewer frames were used when the telescope was observing at low airmass, since 
pupil ghosts produce more severe artifacts when the pupil plane rotates more quickly. These 
ghosts remained in a few of the final science images, but do not affect the central part of 
the array where the QSO is located. Flat fielding was done using Soil's 'special dome flats' 
and custom written reduction script available from the NTT web pages. Registration of the 
dithered images and median co-addition were done with the PHIIRS routines irshift and 
ircoadd, scaling the images to the mode of the sequence. Bright objects were identified in 
the co-added frame and masked out before running through the co-addition process a second 
time. 

The photometric solution was determined in a manner similar to the optical observations 
using IR photometric standards from Persson et al. (1998). Although the March 2004 run 
was photometric, cirrus was present at the end of the October 2003 run. For the affected 
SofI data, we determined K s magnitudes by boot-strapping the photometric solution from 
2MASS point sources. Typically five 2MASS point sources were available per field for the 
bootstrap calibration, the larger errors for the magnitudes of these QSOs reflect the RMS in 
the fit to 2MASS magnitudes 3 A filter transformation 4 was applied in order to place the Soft 



3 One QSO, B0537-286, was observed during both the October 2003 and March 2004. The magnitude 
determined from the bootstrap analysis, K = 16.11 ± 0.04 agrees (within the errors) with the magnitude 
determined under photometric conditions. 

4 http://www. Is. eso.org/lasilla/sciops/ntt/sofi/setup/Zero_Point. html 
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magnitudes on the Persson K s system. However, in the absence of contemporaneous J-band 
magnitudes we have estimated the J— K color as a function of B— K using the dataset of 
FWWOO. Given the scatter in this color-color relation, we expect that the error in assuming 
a J— K color will affect the filter transformation at < 0.02 magnitudes. The final quoted 
i^-band errors in Table 2 include the uncertainty in the fit of the photometric solution, 
photon statistics and the estimated error in the filter transformation added in quadrature. 

Finally, both optical and IR colors were corrected for Galactic extinction according the 
the maps of Schlegel, Finkbeiner & Davis (1998). Our final optical and IR magnitudes are 
presented in Table 2. 

2.3. Literature Magnitudes 

We have supplemented our measurements in Table 2 with magnitudes determined by 
FWWOO, who also achieved near-simultaneous optical and IR (< 6 nights separation) ob- 
servations of a number of CORALS QSOs. In Table 2 we have corrected the K n band 
magnitudes (deduced from the SAAO system) given in FWWOO to K s magnitudes by using 
the conversion determined by 2MASS 5 . We also corrected the photometry of FWWOO for 
Galactic extinction according to the maps of Schlegel et al (1998). 

With these supplemental data, we have optical-IR colors for 45/66 program QSOs, with 
B— K colors for 42. Of these 42, 14 have intervening DLAs; only 3 QSOs with four intervening 
DLAs were not observed by us or FWWOO 6 . The observed sub-sample of CORALS QSOs are 
representative of the whole sample; only scheduling restrictions prevented the observation of 
the full sample. In the remainder of this paper we refer the colors determined from our B 
and K s magnitudes as B— K. 



5 http://www. astro. caltech.edu/ jmc/2mass/v3 /transformations/ 

6 Ellison ct al. (2001) reported 19 intervening DLAs. However, one of these, the z a b s = 1.875 DLA towards 
B2314— 409 has been shown by Ellison & Lopez (2001) to have an N(H I) slightly below the canonical DLA 
limit of 2 x 10 20 cm" 2 . 
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3. Results 

3.1. Comparison of Optical and IR Colors 

The observed B— K color of a given QSO will change with redshift as the intrinsic 
spectrum slides through the filter set. Strong emission features such as Lya and C IV 
can increase the observed flux in a given band, although flux depression due to the Lya 
forest is the dominant source of color shift (e.g. Richards et al. 2003). In Figure 5 we 
demonstrate this effect by calculating the change in B— K as a function of redshift using 
the SDSS composite spectrum of Vanden Berk (2001). The color is calculated relative to a 
reference redshift, in this case z = 2.3. We therefore determine a normalized B— K, (B— K) n , 
which is calculated in the following way. First, we determine the median B— K=2.92 for 
QSOs in the range 2.2 < z < 2.4 since the color changes very little over this interval (see 
Figure 5). Next, we determine a redshift dependent color correction based on the SDSS 
composite spectrum, displayed graphically in Figure 5. In the case of QSOs with intervening 
DLAs, we make a correction (A(B — K) DLA ) for the suppression of flux in the B— band due 
to the Lya line. For each DLA, we simulate a spectrum with a damped profile at the redshift 
z a bs and with iV(H I) given in Table 3. The spectrum is multiplied by the filter transimission 
curve and the resulting magnitude compared with an identical spectrum without a DLA 
absorber. Typical corrections are a few hundredths to a tenth of a magnitude, with a 
maximum A(B — K) DL a = 0.18. The final, normalized B— K color of each QSO is calculated 
by subtracting the median, the redshift dependent term and (in the case of an intervening 
DLA) the absorption correction term [(B - K) n = (B - K) - 2.92 - f(z) - A(B - K) DLA \. 

Table 3 presents the normalized B— K colors for the procedure described above and 
Figure 6 7 demonstrates that the resulting color distribution does not exhibit any redshift 
dependence. There is, however, a large spread in (B— K) n at any given redshift, regardless 
of whether or not a DLA is present, as expected for RLQs (e.g. FWW00). 

3.2. Limits on Extinction 

In order to determine quantatatively whether the colors of QSOs with and without 
DLAs are drawn from the same distribution, we execute a Kolmogorov-Smirnov (KS) test 
on two sub-samples: the 14 QSOs with DLAs and 25 without. The KS test yields a 94% 



7 In this Figure, Table 3 and in the following statistics, we only include QSOs whose B and K s band 
magnitudes have errors less than 0.3 magnitudes. In practice, this excludes only 3 of the FWW00 QSOs 
which have poorer photometry. Our conclusions are unaffected by the exclusion of these three. 
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probability that the two sub-samples are drawn from the same distribution. The mean 
(median) normalized colors of QSOs with DLAs are +0.12 (+0.01) and —0.10 (+0.02) for 
QSOs without DLAs with RMS scatters of 0.56 and 0.47 respectively. A student's T-test 
yields a test statistic = 1.32, indicating that the mean colors of the two samples are not 
significantly different. 

In order to quantify a limit to the amount of reddening associated with the DLAs in 
our sample, we run a series of simulations that artificially include/remove reddening in the 
spectra of QSOs with DLAs. We then calculate the KS statistic compared with the no DLA 
sample. In each simulation, we redden the observed B— K color of those QSOs with DLAs 
assuming a fixed EiB — V) for each DLA for both Galactic and SMC extinction curves. We 
use the parameterizations of SMC (R v = 2.93) and Galactic (R v = 3.08) extinction from 
Pei (1992) and a second functional fit to Galactic reddening from CCM89, as described in 
section 1.3. The colors of QSOs without DLAs remain unchanged in these simulations. 

To demonstrate how a small change in E(B — V) can affect the color distribution and 
KS statistic we show the B— K colors for a range of simulations in Figures 7 and 8. The 
corresponding KS probability is given above each panel. The results of the full simulation 
gamut for the three extinction curves are shown in Figure 9, where a positive E[B — V) 
indicates that the spectrum was reddened, whereas a negative value indicates that the effects 
of dust were removed. These simulations rule out E(B — V)> 0.04 (SMC) and E(B — V)> 
0.07 (MW) at 99% confidence (i.e. KS probability < 0.01 for E(B - V) < -0.04, -0.07, 
with no difference between the Pei and CCM89 Galactic models). 

Alternatively, we can consider a reddening that depends on the iV(H I) of each DLA. 
Reddening-to-gas 8 ratios have been determined in local galaxies, e.g. 

f(b v) N{HI) m 

E{B ~ V) - 4.0 x 10^ (1) 
for the SMC (Bouchet et al. 1985) and 

E(B-V)= N{HI l (2) 
v ; 5.8 x 10 21 V ; 

for the Milky Way (Bohlin, Savage & Drake 1978). 



8 Some authors refer to the relationship between E(B — V) and 7V(H I) as a dust-to-gas ratio. However, 
we prefer to use the term reddcning-to-gas ratio to avoid confusion with the parameters k and k that are 
often used in the literature to denote dust-to-gas ratios inferred, for example, from chemical abundances. 
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If reddening is removed from the DLA-QSO sample in accordance with the SMC reddening- 
to-gas relation in Equation 1 the KS probability is 0.39, an inconclusive result. Correcting 
the colors of QSOs with DLAs using the more extreme Galactic relation in Equation 2 results 
in colors that are significantly bluer than the no-DLA QSOs, i.e. the reddening correction 
is too large. The KS probability in this case is 0.03, ruling out a Galactic reddening-to-gas 
relation at 2 a. 

In reality, reddening is likely to be more stochastic than the simple dust recipes adopted 
here, and extinction curves are likely to vary from DLA to DLA. Our dataset is currently too 
small to quantify the reddening-to-gas ratios for DLAs through modelling techniques such 
as x 2 minimization, although we can rule out a Galactic scaling relation at 97% confidence. 
However, future work with larger samples, e.g. with DLAs identified from the SDSS, will 
be able to investigate the reddening-to-gas relation in more detail. Finally, we note that a 
detailed treatment of other such effects as non-negligible intergalactic medium reddening, 
host galaxy contributions or exotic extinction laws are beyond the scope of this paper. 

3.3. X-ray Observations of the Reddest QSOs 

To further investigate the 3 QSOs with the reddest B— K colours, we have searched for 
archival X-ray observations of these sources. Given a sufficient signal-to-noise in the X-ray 
data, it is possible to investigate whether 'intrinsic' absorption (i.e., associated with the 
QSO) is present. In principle the X-ray data can measure, or place constraints upon, the 
integrated hydrogen column along the line of sight (and by assuming a dust-to-gas ratio, on 
the reddening), as well as the redshift at which the absorption takes place. Since the inferred 
iV(n I) column densities of intrinsic absorbers are frequently an order of magnitude higher 
than typical intervening DLAs, it is plausible that significant reddening may be induced by 
such an absorber (although an important caveat is the large possible range in the reddening- 
to-dust ratios). 

Of the three sources with (B— K)„ > 1 we found that X-ray observations had been 
obtained for the QSOs B0438-436 and B0458-020, both of which exhbit DLA absorption. 
Unfortunately, no X-ray observations have been obtained for B1318— 263, the system without 
a DLA absorber. 

B0438— 436 has long been known to exhibit significant excess (i.e. above the Galactic 
contribution) absorption. Evidence for this excess absorption has come from the analysis of 
ROSAT, ASCA and XMM observations (Serlemitsos et al. 1994; Elvis et al. 1994; Cappi 
et al. 1997; Brocksopp et al 2004). The best quality data so far are those obtained with 
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XMM and the best-fit model (Brocksopp et al 2004) is consistent with a solar metallicity 
absorber with N(H) = 1.3 x 10 22 cm -2 , which Brocksopp et al. assumed to be at the redshift 
of the QSO 9 . If we assume that ionized and molecular gas is negligible, this column density 
translates into a E(B — K) ~ 18.5 and ~ 5, for a Galactic and SMC reddening-to-dust ratio 
respectively and extinction curves from Pei (1992). These values are considerably higher 
than the E(B — K) observed for this source (see Figure 6). While there are a number of 
explanations for this discrepancy (e.g. dustless gas near the QSO or AGN extinction curves 
that are significantly different from Local Group galaxies, e.g. Gaskell et al. 2004) a large 
reddening due to assoaciated absorption is very plausible. 

We have retrieved archival Chandra ACIS data for B0458-020 (ObsID 2985) with a 
total exposure of 77.76 ksec. The data were analysed in the 0.3-10.0 KeV energy range using 
CIAO 3.2.1 software and the latest calibration files. B0458— 020 suffers from relatively high 
Galactic absorption, corresponding to N(H I) = 7.5 x 10 20 cm~ 2 , which hinders the detection 
of intrinsic absorption which is redshifted to the soft (low energy) part of the spectrum. 
Indeed, a spectral fit assuming a simple power-law spectrum with only Galactic absorption 
yields a good fit to the data, i.e. Xmd = 0-8 for 155 degrees of freedom. The best-fit photon 
index V (/„ oc z/ r+1 ergs s" 1 cm -2 ) was found to be 1.55±0.02 (see Figure 10). The inclusion 
of absorption at either the QSO or the DLA redshift gave no significant improvement to the 
fit: N(H) < 2 x 10 21 cm" 2 with 99% confidence (F-test statistics F ~ 0.002 for A^ 2 ~ 7.3 
and 1 additional parameter). Using the Pei (1992) extinction laws this corresponds to an 
upper limit E(B — K) < 2.8 and < 0.7, for a Galactic and SMC reddening-to-dust ratio 
if the absorber is located in the host galaxy. Associated absorption could therefore be the 
cause of the red color of this QSO, although better X-ray data would be required to confirm 
this. However, we note that Carilli et al. (1998) found that 80% of red QSOs in their sample 
exhibited strong 21 cm absorption at the QSO redshift, supporting the view that red colors 
may be often caused by intrinsic absorption. 



4. Discussion 

In a series of papers, Fall & Pei (1989), Fall, Pei & McMahon (1989) and Pei, Fall & 
Bechtold (1991) investigated whether the spectral indices of QSOs with DLAs were system- 
atically steeper that QSOs with no intervening galaxy. In particular, in the latter of these 



9 With sufficiently high quality data, it is possible to determine the absorption redshift. However, asso- 
ciated absorption is often the working assumption since the column densities greatly exceed that of typical 
intervening absorbers. In this case, the intervening DLAs has 7V(H I) = 6 x 10 20 cm -2 . 
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papers, new data with higher spectrophotometric accuracy were obtained and a sample of 
66 QSOs analysed for reddening. The distribution of spectral indices in the DLA sample 
(20 QSOs) was found to differ at the 99.999% confidence level from the no-DLA (46 QSOs) 
sample. The inferred dust:gas ratio was 5-20% of that in the Milky Way. Further support 
for observed reddening was reported by Outram et al. (2001) who found redder B— R colors 
for QSOs with 0.5 < z < 1.5 high equivalent width Mg II absorbers. The inferred color 
excess was E(B — V) ~ 0.04. More recently Wild & Hewett (2005) have found a convincing 
reddening signature associated with low z Ca II absorbers in the SDSS, although still at 
a modest level: E(B — V) = 0.06. Individual cases of DLA reddening have been reported 
based either on differential reddening in lensed systems (e.g. Zuo et al. 1997 for a DLA at 
Zabs = 1-4) or the detection of the 2175 A graphite feature (Junkkarinen et al. 2004 for a 
DLA at z a bs = 0.5). However, some caution is required when comparing inferences on dust 
extinction at different redshift ranges since the amount, source and composition of dust may 
differ, for example as the dominant source of dust shifts from Type II supernovae to lower 
mass stars (e.g. Hirashita et al. 2005). Nonetheless, these works apparently provide robust 
evidence for the systematic reddening of QSOs with intervening absorption systems at both 
high and low redshifts. 

Recently, however, Murphy & Liske (2004) re-visited the spectral fitting approach with 
a large statistical sample drawn from the SDSS. In addition to its size (~ 1450 QSOs, with 
72 DLAs) this sample has the advantage that its data acquisition is homogeneous and the 
spectral indices are fit over a wider wavelength range than was done in earlier work. With 
these improvements, Murphy & Liske (2004) found no steepening of the spectral indices for 
QSOs with DLAs and determined a limit of E(B - V) < 0.02 (99% confidence). Although 
Murphy & Liske adopt a similar approach to Fall, Pei and collaborators in fitting spectral 
indices, they do so with an order of magnitude larger sample, wider wavelength coverage 
and uniform data acquisition. Since QSOs have a wide range of intrinsic colors and spectral 
shapes, using the broadest baseline possible is important for identifying a systematic dust 
reddening. 

For the first time, optical-IR colors can now contribute to the on-going debate. Although 
our sample is modest in size in comparison with the SDSS (although comparable with the 
samples of Pei, Fall and collaborators), the broad wavelength coverage provides significant 
leverage for detecting reddening. For example, whereas the B— K color changes by 0.74 for 
an E(B — V)=0.05 (SMC extinction z cm = 3.0, z^ s = 2.6), the B — R color changes only by 
0.35. Moreover, the sample of QSOs on which we conduct this work is based on an optically 
complete survey of radio-loud quasars. Therefore, our sample of DLAs should not be subject 
to selection effects associated with dust extinction. 
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Our results are presented graphically in Figure 6 where we show the normalized B— K 
colors for QSOs with and without DLAs. The mean color of QSOs without DLAs is —0.10, 
compared with +0.12 for the DLA sub-sample. Both color distributions have RMS scatters 
~ 0.5. Neither a student's T-test nor a KS test find any significant difference in the B— K 
colors of the two sub-samples. We determine a 3<r limit of E(B — V) < 0.04, 0.07 for 
SMC and Galactic extinction respectively based on the results of our KS test simulations. 
Although not quite as sensitive as limits from the SDSS (Murphy & Liske 2004), our sample 
is optically complete and could potentially contain heavily extinguished QSOs which would 
not have been fond in magnitude limited samples. However, the low E(B — V) allowed by 
this work shows that even in optically complete samples, the amount of reddening is small. 
Such low values of EiB — V) are consistent with the small H2 fractions observed in DLAs; 
whilst a handful of absorbers have log f(H 2 ) ~ —2, most have log f(H 2 ) < —5 (Ledoux, 
Petitjean & Srianand 2003). Indeed, Tumlinson et al. (2002) have shown that low molecular 
fractions are common in Galactic and Magellanic Cloud sightlines when E(B — V) < 0.08. 

Although this work has not found any indication that QSOs are significantly reddened 
by intervening DLAs, evidence that some dust is present in DLAs is well established from 
the abundances of refractory versus volatile elements such as Cr and Zn (e.g. Pettini et al. 
1994, 1997). Akerman et al. (2005) have obtained Zn and Fe abundances for many of the 
DLAs in our sample; we reproduce these abundance results in Table 3 adopting the solar 
reference values from Lodders (2003). Vladilo & Peroux (2005) have recently suggested that 
reddening will increase linearly with N(Zn) and become especially severe when log N(Zn) 
> 13.2. In Figure 11 we plot the normalized B— K color versus N(Zn) using the column 
densities reported in Akerman et al. (2005) 10 . We detect no trend between normalized color 
and log N(Zn), although we have only 3 DLAs with log N(Zn) > 13.0. Nonetheless, Akerman 
et al (2005) find that the mean weighted metallicity of CORALS DLAs is only marginally 
higher (although consistent within the error bars) than previous samples. 

In Figure 12 we investigate whether there is any correlation between N{Yi I) or deple- 
tion (as measured by [Fe/Zn]) and reddening. The two reddest QSOs which also show DLA 
absorption have slightly larger than average ([Fe/Zn] = —0.3, Prochaska & Wolfe 2002), but 
still modest depletions. The rest of the sample (including two systems with comparable de- 
pletions to the reddest DLA-harboring QSOs) do not show any correlation between depletion 
and reddening. Similarly, the N(HI) of one of the three very red QSOs in our sample has the 



10 We make the usual assumption that Zn II is the dominant ionization stage of this element and that the 
contribution from other ionzization stages is negligible. Moreover, in Figures 11 and 12 we include 2 DLAs 
within 3000 km s _1 of the QSO, plus one further absorber whose N(R I) is just below the canonical DLA 
limit. 
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DLA with the highest iV(H I) amongst CORALS absorbers. However, a second high iV(H I) 
DLA has apparently not reddened its background QSO to a similar extent and the reddest 
QSO with a DLA has only a modest iV(H I). Overall, there is no obvious trend between 
(B— K) n and iV(H I). On the other hand, Khare et al. (2004) do report tentative evidence 
for a correlation between the inferred E(B — V) (based on broad band SDSS colors) and 
[Cr/Zn] (analogous to the [Fe/Zn] plotted in Figure 12). However, the Khare et al. result 
is based on the two systems in their sample with [Cr/Zn] < 1, and only one of the DLAs in 
our sample may have such a high dust:gas ratio. The normalized B— K color of this QSO 
is mildly redder than the median color of the other QSOs, although this can certainly not 
be taken as evidence for a correlation between depletion and reddening. It would be of 
great interest to further investigate the correlation between reddening and dust:gas ratio in 
a larger sample of objects. 

Prantzos & Boissier (2000) have suggested that dust bias is responsible for the observed 
anti-correlation between [Zn/H] and A(H I), and propose an empirical 'dust filter', i.e. that 
DLAs in optically selected samples all have log iV(H I) + [Zn/H] < 21. In Figure 13 we 
plot iV(H I) versus [Zn/H] from the compilation of DLAs in Kulkarni et al. (2005); the 
dust filter suggested by Prantzos & Boissier (2000) is shown by the dashed line. We can 
combine the suggested cut-off of log iV(H I) + [Zn/H] < 21 with a reddening-to-gas ratio 
which scales with metallicity in order to calculate what E(B — V) is implied. For example, 
it is straightforward to show that if the Galactic reddening-to-gas ratio given in Eqn 2 is 
scaled by [Zn/H] and satisfies the criterion log JV(H I) + [Zn/H] < 21, then E(B — V) < 
0.17. This is the E(B-V) value required to ensure that no high N(HI), high [Zn/H] systems 
are observed; E(B — V) = 0.17 is indeed considerably larger than the limits determined by 
this work and that of Murphy & Liske (2004). We can alternatively calculate where the 
observed cut-off in Figure 13 would lie for a given assumed E(B — V). For example, an 
E(B — V) < 0.05 implies log iV(H I) + [Zn/H] < 20.46 (again using Equation 2 scaled by 
i(][ Zn / H l); this cut-off is shown in Figure 13 by a dot-dash line. If future work determines 
an E(B — V) « 0.05, it may (depending on the assumed reddening-to-gas ratio) become 
difficult to quantitatively reconcile the N(R I)-[Zn/H] distribution with dust bias. 

In Figure 13 we also plot the iV(H I) and [Zn/H] values for the CORALS DLAs, which 
should not be subject to dust bias. The figure shows that CORALS DLAs are not distinguish- 
able from other (literature) DLAs in terms of their N(R I)-[Zn/H] distribution. Moreover, 
they do not lie above the empirical cut-off which has been proposed to be due to dust. How- 
ever, in the chemical evolution models of Churches et al.(2004), < 9% of DLAs (depending 
on the choice of spin parameter) are expected to be in this 'dust forbidden' region (A. Nelson, 
2005, private communication). It therefore seems premature at this point to interpret the 
lack of high iV(H I), metal-rich CORALS DLAs as evidence that the anti-correlation is not 
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caused by dust. 

If future work concludes that dust is the reason for the anti-correlation in Figure 13, then 
the locus of data points compared with the theoretical lines and the lack of high CORALS 
values may indicate that typical E(B — V) values are simply intriniscally low in most DLAs. 
Indeed, Ellison, Kewley & Mallen-Ornelas (2005) found E(B - V) < 0.2 even in the central 
parts of absorber counterparts at z < 0.5. Alternatively, the implied E(B — V) values from 
this and the work of Murphy & Liske (2004) may be explained if the dust extinction is 
relatively grey, as been proposed to explain the E(B — V) values observed towards some 
gamma-ray burst hosts (e.g. Vreeswijk et al. 2004). 

One of the greatest uncertainties in this work is the unknown nature (i.e. composi- 
tion, grain size and, therefore, extinction curve) of dust in damped systems. The main 
observational distinction between SMC and Galactic extinction curves (with the LMC being 
somewhat intermediate between the two, see Figure 1) is the presence of a 'bump' at 2175 
A, attributed to graphite. A handful of detections of this feature at cosmological redshifts 
exist associated with 3 Mg II absorbers with z ahs ~ 1.5 in the SDSS (Wang et al. 2004), a 
high N(HI) DLA at ^ abs ~ 0.5 (Junkkarinen et al. 2004) and several cases of galaxy lenses 
(e.g. Falco et al. 1999; Toft et al. 2000; Motta et al. 2002; Wucknitz et al. 2003; Munoz et 
al. 2004). However, gravitationally lensed systems are clearly special cases since the impact 
parameter through the galaxy is necessarily small. Malhotra (1997) have claimed to see the 
2175 A feature in stacked (non-lensed) QSO spectra exhibiting Mg II absorption, although 
stacks of SDSS spectra have failed to reproduce this (M. Murphy, 2005, private communica- 
tion; Menard 2005; Khare et al. 2005; Wild & Hewett 2005). Although further investigation 
is warranted, this may be evidence that Galactic dust is not generally applicable for absorbers 
towards unlensed QSOs. Theoretical models of dust in high redshift galaxies enriched by 
high mass Type II SN may be more relevant than local extinction curves. However, the 
extinction properties amongst these models also vary tremendously depending on progenitor 
mass and mixing (e.g. Hirashita et al. 2005). 



5. Conclusions 

In this work we have investigated whether dust in DLAs can cause a significant, sys- 
tematic reddening of background QSOs by obtaining B— K colors for an optically complete, 
radio-selected quasar sample. We determine normalized colors which account for redshift de- 
pendence (analogous to a K-correction which accounts for Lya flux decrement and emission 
lines) and suppressed B-band flux due to DLA absorption. The mean normalized B— K color 
of the DLA sub-sample is 0.12±0.56, compared with — 0.10±0.47 for the no-DLA sub-sample. 



-17- 



Both a student's T-test and KS test indicate that there is no statistical difference between 
the distribution of B— K colors in the two sub-samples. We place a limit on the amount of 
reddening that may be present, by correcting the colors of QSOs with DLAs using a range of 
E(B — V) values. Adopting an SMC extinction curve, we place a 3a limit on the systematic 
reddening: E(B — V)< 0.04, assuming a fixed reddening for each DLA. In general, there is 
no correlation between depletion, N(Zn) or iV(H I) and reddening. Finally, we have searched 
for archival X-ray data of the three reddest QSOs in our sample in order to test whether 
their extreme (B— K)„ colors maybe due to intrinsic (i.e. at z ~ z cm ) absorption. Of these 
three QSOs, two have data available, and the spectrum of one (B0438— 436) shows evidence 
for a high N(R I) absorber associated with the QSO which may be the cause of the very red 
color. 

These results support relatively low extinctions towards high redshift DLAs, and agree 
with the E(B — V) values inferred from the optically selected SDSS work of Murphy & 
Liske (2004). Therefore, although a handful of individual absorbers may cause more severe 
reddening (e.g. Junkkarinen et al. 2004), these seem to be the exception rather than the 
rule. A number of observational works have also inferred relatively small reddening at lower 
redshifts (e.g. Khare et al. 2005; Menard 2005; Wild & Hewett 2005). Combined with 
previous work which has found that the number density and gas content of intervening 
absorbers (Ellison et al. 2001, 2004) and their metallicities (Akerman et al. 2005) are 
not largely different from samples whose optical completeness exceeds B ~ 19, it seems 
increasingly unlikely that dust bias has severely skewed our view of high redshift DLAs. Put 
another way, the red(dened) QSO population is unlikely to be linked to dust in high redshift 
DLAs. 
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Table 1 — Continued 



QSO 


4m 


Observing 


Integration 


Integration 


SuSI2 






Run 


Time Soil (s) 


Time SuSI2 (s) 


Filter 
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Table 2. CORALS QSO magnitudes 
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16 


71±0 


04 


This 


work 


B1010- 


-427 


2 


954 






17 


25±0.07 


16.61±0.03 
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Table 2 — Continued 
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z 
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Table 2 — Continued 



QSO z em Bj UKST B V K Ref 



B2314-409 2.448 19.0 18.26±0.06 18.03±0.05 15.58±0.08 FWW00 

B2315-172 2.462 21.0 

B2325-150 2.465 20.0 

B2351-154 2.665 18.7 18.82±0.07 18.60±0.06 16.37±0.26 FWW00 



Table 3. Normalized Colors and Intervening Absorber Properties for CORALS QSOs with 

Near- Simultaneous B— K Colors 



QSO 


z 


em 
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Table 3 — Continued 



QSO 


4m 




N(HI) 


(B - K) n 


[Zn/H] 


[Fe/H] 


Abundance Ref. 
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-0.32±0.06 
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Fig. 1. — Extinction curves for the parameterizations of Pei (1992) and Cardelli, Clayton & 
Mathis (1989). Assumed R v values are 2.93 (SMC), 3.16 (LMC) and 3.08 (MW). There were 
no data available at wavelengths < 1000 A, hence the divergence in the FUV for the Galactic 
curves represents an extrapolation of the fitting functions. Filter transmission curves along 
the base of the plot illustrate what rest wavelengths correspond to the B and K s bands for 
z = (solid) and z = 3 (dashed). 
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E(B-V) 



Fig. 2. — Theoretical reddening of a z em = 3.0 QSO with an intervening absorber at z a hs = 
2.0, 2.3, 2.6, 2.9. We assume a Galactic extinction curve based on the parameterization of 
Pei (1992). 
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E(B-V) 

Fig. 3. — Theoretical reddening of a z em = 3.0 QSO with an intervening absorber at z^s = 
2.0, 2.3, 2.6, 2.9. We assume a Galactic extinction curve based on the parameterization of 
Cardelli, Clayton & Mathis (1989) with R v = 3.08. 
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E(B-V) 

Fig. 4. — Theoretical reddening of a z em = 3.0 QSO with an intervening absorber at z^s = 
2.0, 2.3, 2.6, 2.9. We assume an SMC extinction curve based on the parameterization of Pei 
(1992). 
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Fig. 5. — Change in B— K color as a function of redshift determined using the SDSS composite 
spectrum of Vanden Berk (2001). Color changes are relative to z = 2.3. 
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Redshift 

Fig. 6. — Redshift distribution of normalized B— K colors for QSOs with (filled circles) and 
without (open stars) DLAs. 
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E(B-V) = 0.01, KS probability = 0.251 E(B-V) = 0.02, KS probability = 0.084 




2 2.5 3 3.5 4 2 2.5 3 3.5 

Redshift Redshilt 

E(B-V) = 0.03, KS probability = 0.013 E(B-V) = 0.04, KS probability = 0.000 
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Redshift Redshift 

Fig. 7. — Redshift distribution of normalized B— K colors for QSOs with (filled circles) and 
without (open stars) DLAs with reddening added to the former population assuming a Pei 
(1992) SMC extinction law. Each of the four panels illustrates a different degree of reddening 
and shows the KS probability that the DLA and non-DLA QSOs are drawn from the same 
color distribution. 



- 36 



E(B-V) = -0.01, KS probability = 0.995 E(B-V) = -0.02, KS probability = 0.234 




2 2.5 3 3.5 4 2 2.5 3 3.5 

Redshift Redshilt 

E(B-V) = -0.03, KS probability = 0.047 E(B-V) = -0.04, KS probability = 0.011 
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Redshift Redshift 

Fig. 8. — Redshift distribution of normalized B— K colors for QSOs with (filled circles) and 
without (open stars) DLAs with reddening removed to the former population assuming a Pei 
(1992) SMC extinction law. Each of the four panels illustrates a different degree of reddening 
and shows the KS probability that the DLA and non-DLA QSOs are drawn from the same 
color distribution. 
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Fig. 9. — KS probability that the (B— K) n colors of QSOs with and without DLAs are 
drawn from the same distribution, for simulations which include reddening in the former 
population. Positive values of E(B — V) indicate that the effect of dust (i.e. reddening) 
was added to the simulated spectra and negative values indicate that the effect of dust was 
removed. The vertical dotted line indicates that the raw (B— K) n have been used, i.e. no 
reddening included. 
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Fig. 10. — Best-fit y 2 plots of the B0458— 020 Chandra observations. The fit and sigma 
((data - model) / error) plots are shown after convolving with the instrument response 
(top and bottom panels, respectively). The data were fitted with a single power-law model 
absorbed by a Galactic foreground iV(H I) column (thin continous line in top panel). 
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Fig. 11. — Normalized B— K color of CORALS DLAs as a function of Zn column density. 
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Fig. 12. — Top panel: Normalized B— K color as a function of depletion, as measured by 
[Fe/Zn]. Bottom panel: Normalized B— K color as a function of neutral hydrogen column 
density, N(HI). 
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Fig. 13. — Metallicity versus iV(H I) column density for DLAs taken from the compilation 
of Kulkarni et al. (2005) and for the CORALS DLAs (Akerman et al. 2005), shown by open 
and filled circles respectively. The dashed line shows the proposed 'dust filter' of Prantzos 
& Boissier (2000): log JV(H I) + [Zn/H] < 21, which corresponds to E(B — V) < 0.17 for 
a Galactic gas-to-reddening law scaled by metallicity. The dot-dashed line shows where the 
cut-off would lie if E(B - V) < 0.05, i.e. log JV(H I) + [Zn/H] < 20.46. 



